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Abstract Recently a formalism for a direct treatment of the Faddeev equa- 
tion for the three-nucleon (3N) bound state in three dimensions has been 
proposed. It relies on an operator representation of the Faddeev component 
in the momentum space and leads to a finite set of coupled equations for 
scalar functions which depend only on three variables. In this paper we pro- 
vide further elements of this formalism and show the first numerical results 
for chiral NNLO nuclear forces. 

"o ' 

5 ' 1 Introduction 

Chiral three-nucleon (3N) forces at next-to-next-to-next-to- leading order (N3LO) 
(see for example Refs. [IJ[2]) have a very rich spin and momentum structure, 
so calculations of the 3N system based on standard partial wave representa- 
tions are very tedious [3] . Recently, the Faddeev equations for the 3N system 
have been formulated directly with vector variables [H[MM7]. in contrast to 
Refs. [HG3|6], where spin and isospin couplings are treated explicitly, the key 
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ingredient in Rcf. [7] is an operator representation of the 3N bound state. 
Such an operator form was given long time ago in Ref. [8 and later re-derived 
in Ref. [TO]. 

Following this formulation, one represents two-nucleon (2N) t-operators, 
3N forces and Faddeev components of the 3N wave function as linear com- 
binations involving products of scalar functions and scalar operators built 
from the spin operators and momentum vectors. By simple manipulations it 
is later possible to remove the spin operators, which leads to relations among 
scalar functions of momentum vectors only. In this paper we provide further 
elements of this formalism. 

An alternative formulation of the 3N bound state equation is discussed in 
Sect. 2. There we introduce an additional set of scalar coefficients required 
by the modified integral kernel. 

In Sect. 3 we consider the full 3N bound state and show how its scalar 
coefficients can be obtained from the scalar functions used to represent the 
Faddeev component. In the same section we discuss the normalization of the 
full wave function, as well as the expectation values of the kinetic energy and 
the 3N potential energy. 

Details of our numerical performance and some numerical results are pro- 
vided in Sect. 4, where we use a simpler and more efficient way of introducing 
3N forces. We end with a brief summary and outlook in Sect. 5. 

2 An alternative form of the Faddeev equation using the operator 
expansion 

For the convenience of the reader we first remind of the building blocks of 
our formalism. 

Working with the isospin concept, nucleons are treated as identical parti- 
cles so one needs just one equation for the Faddeev component \tp). In Ref. [7] 
this equation is written in the following form 

IV) = G tP\^j) + (1 + G Q t)G Q V^(l + P)\i>) , (1) 

where Go = e~h 1S * ne ^ ree 3N propagator for the internal 3N energy 
E, t is the 2N transition operator for nucleons in the (2, 3) subsystem and 
P = P12P23 +P13P23 is the permutation operator. The term is that part 
of the 3N force, which is symmetrical under exchange of nucleons 2 and 3. 
Of course, the choice of the subsystem pair (2, 3) is arbitrary. 

Alternatively, the Faddeev equation can be written without the transition 
matrix element t as 

IV) = G V(1 + P)\i>) + G y (1) (l + P)\i>) , (2) 

using directly the 2N potential V for nucleons 2 and 3. 

In order to account for charge independence and charge symmetry break- 
ing, matrix elements of the 2N potential between the eigenstates of the 2N 
isospin operator are given as 

(tm t I V I tW) = Sw <W-v ytmt ■ ( 3 ) 
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It is well known [TT] that the most general 2N force V tmt has in the momen- 
tum space the following form 



(P I v tmt I P') 



Etrn t 

i=i 



(p,p') w,-(er (2) ,<T (3) ,p,p'), 



(4) 



with v- mt (p, p') being scalar functions of the initial (p') and final (p) rela- 
tive Jacobi momenta. We choose the same spin-momentum scalar operators 
Wj(<T( 2 ), C(3), P, p') as in Ref. [T2] 



Wi(er (2) ,er (3) ,p, p 
U? 2 (t7-(2),0-(3),p ! P' 
W3(o-(2),<T (3) ,p, p' 

W4,((T( 2) ,cr {3) ,p, p 
w 5 (er(2),o-(3),P, p' 
^6(0-(2),0-( 3 ),p, p 



) = 1, 

) = f (2) • <T(3) , 

) = -i(o-(2) +cr(3) ) • (P x P') , 

) = CT(2) ' (P X P')CT(3) ' (P X P') 

) = <T(2) ■ (p + p')^) • (p + p') 

) = CT(2) ' (P - P') Cr (3) ' (P - P') 



(5) 



(Note this set is different from the one used in [7].) Also for the transition 
operator t 1 , generated by the Lippmann-Schwinger equation, 



Lrrit 



ytm t + ytm t Q^ t 



trn± /^i2N j-trrit 



(6) 



the same expansion is valid but the corresponding scalar functions i* m *(p, p') 
depend also on the 2N energy, E 2 n- 

In the next step, Up), V, t and are expanded into the three possible 
3N isospin states, | (i^) T), given in terms of the subsystem isospin t and 
the total 3N isospin T: 



V 



= E 

tT 

E 

tT,t'T' 



*= E 



tT,t'T' 

E 

tT,t'T' 



4 ) 7 )!''"' 



t-\T)V t 



(1) 



tTA'T 




Consistent with ([3]), we assume that 

VtT,t'T' — $tt> VtTT' 



(7) 
(8) 
(9) 
(10) 

(11) 



ttT,t'T' — Stt' ttTT' 



(12) 
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We consider only such 3N forces which conserve the total 3N isospin: 

V&] VT ,=&rr>V$> T . (13) 
The isospin matrix elements of the permutation operator P read [13] 

( (4) T 1 p 1 (4) r ' } = 5 TT'F t MP^ n P% 1 + (-)* + *>iT^T) . (14) 

where Ftt'T are geometrical factors and Pf™ acts only in the spin and momen- 
tum spaces. Using the above information, the Faddeev equation ([2]) becomes 

\M = g J2 v *tt' J2 l 5 "' + F "' T ' {p^PiT + i-) t+t ' PlTPif)) \M 

T' V 

+ GoJ2 v $t E (<**'*" + Pt'fT (p™Pk + (-) 4 ' +t '>iT^T)) 

(The isospin projected equation ^ can be found in Ref. 7 ].) As expected, 
the charge independence and charge symmetry breaking in the 2N force leads 
to the coupling of total isospin T = \ and §. 

In the following we will consider the above equation in the momentum 
space spanned by the two relative Jacobi momenta p and q. In Eq. ([T]) we 
encounter the off-shell 2N i-matrix 

(pq|ttT7* Ip'q') = <Hq - q%TT' (p, p', E q ) (16) 

with E q = E 2 n = E — ^g 2 - In the case of Eq. @ we only need the 2N 
potential matrix elements 

(pq\V tT T> |p'q') - 5(q - 4)V t TT> (p, p') . (17) 
The matrices for the permutation operators read |14j 

(pqlAT^TIpV) = S(p - 7r(q, q'))*(p' - 7r'(q, q'))^^ (18) 

(pq| AT^TIpV) = s(p + *(q, q')W + w'(q, q'))^^ (is) 

with 

7r(q,q')-^q + q' (20) 

7r'(q,q') = -q-^q / . (21) 

As pointed out in Ref. [7] , it is crucial for the present formulation to write 
the Faddeev amplitude iptr(p, q) = (pq | V'tr) in the operator form (TU]. It 
reads 

8 

^T(p,q) = E^T(P>q) C 'ilx m ) I (22) 
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where \x m ) — \{^\)\ m ) is a specific state in which the three spin-1/2 states 
are coupled to the total angular momentum quantum numbers of the 3N 
bound state. The functions </>[^(p, q) = </>[^(|p|, |q|, p-q) are scalar functions 
and the operators Oi are given as 



01 = 1, 

02 = -j= <t(23) ■ <r (1) , 



3 i 



°3 = \/--<T (1) .( P xq) 



4 

o 5 = - 



1 

f 



i<t{23) • (p X q) - (<x (1) x <r(23) ) ■ (p x q) 



<r(23) • (p x q) - - (<r (1) x <r(23)) • (p x q] 



° 6 = vf ( <T ^ 23 ) ' p ' p ~ \ <T ( 23 ) ' ' 
° 8 = § TTg ( <T ^ 23 ) ' q ' p + cr ( 23 ) ' p^m ' 3 



-p • qer 



(23).<t (1) ) 



(23) 



where cr(23) = 5(072) — 0(3) )• We use unit vectors, p and q. Note that we 
actually use the second set of operators from Ref. [7j . Note also that several 
misprints in Eq. (31) of Ref. [7] has been now corrected using [10 . 

As already mentioned, we can use the most general operator structure of 
the 2N t-matrix t t TT' (p, p', E q ), which is given by 



t t TT'(p,p',E q ) = ^t$ T ,(p,p',Eq) w,(<x (2) ,<T (3) ,p,p'). (24) 
i=i 

Here the matrix elements 4tt' (P»P' > Eq) are scalar functions and the oper- 
ators Wj(<T( 2 ), <T(3), p, p') are given in Eq. (J5]). 

In principle, the 3N force operator V^/^(p, q, p', q') also allows for an expan- 
sion 



^tt'HP) q> p'' <*') = w «'t(p> q> p'> q') n i ( CT (i) > CT (2) , 0(3) , p, q, p', q'X 2 5) 



with a strictly finite number of terms [5] but for computational reasons we 
treat it as one object and do not use any expansion. This and the fact that our 
integral kernel is simpler than in Ref. [7] , leads to a substantial simplification 
and a reduction of the number of scalar coefficients introduced in Ref. [7] . 
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Using exactly the same algebra as in Ref. [7], we arrive at the final form 
of the Faddeev equation ©, which corresponds to Eq. (50) in Ref. [7J: 

8 

E (pq> pq) 4>tT (Pi q) = 
i=i 

*-* m 4m y 17 T< k' = l 

6 

E ( u tTT' (P> ^fa ^)) D kok' (pq, P7r(q, q'), 7r'(q, q')q') 

+ (p, -7r(q, q'))- ^' (pq' pf-^fa q'))» w'(q. q')q') 



+ 



+ 



E^'^'^K(q,q')q' 

~ 6 8 

/ d V E E E v tTT' (p, p')^fcifc' (pq^ pp'. p'q)^ (p', q) 

/ dVA'E E <' T (pq,pqp , q / ,p'q , )^T ) (P , ,qO 

^ f fc'=i 

f AWE E (G^(pq,pq^(q',q")q'^'(q',q")q") 



1 


E- 


Pl _ 

m 

l 


4m y 


E- 


m 

l 


4m y 



E _Pi_ _3_„2 ^ 

m 4m " £' fe — 1 

+ (-) t 'G;»' T (pq, pq(-7r(q', q"))q', 7r'(q', q")q" 



x^M»T^'(q'q"),q"). (26) 

We use here seven sets of scalar coefficients: Cfej, Ayi, ^jj, E ki T , 

G k - T and G kk , T . Cm, ^>kjk' and £> fej - fe , are exactly the same as in Ref. [7J, 
where they appear in the GgtP | tp) part of the Faddeev kernel. The co- 
efficients E l kk T , G kk T and G kk , T , which are required to deal with the 3N 
force in the Faddeev equation, are slightly modified versions of the quan- 
tities appearing in Ref. [7J, because we do not expand V^)p(p, q, p', q') = 

v 'i ( t'T(' T (i)> cr (2); cr (3) J p,q, p',q') : 
£££ T (pq»pqp'q',p'q') = 

^2(x m \Ok(cr(i), o-( 2 ), o-( 3 ),p, q) 

m 

xl/ t ( t ^(cr (1) ,cr (2) ,cr (3) ,p,q, p',q') O k > (<X(i), <r (2) , er (3) , p', q')|x' m ) , (27) 
G^(pq,pqp'q',p"q") = 

~~ 2 ( x ™ l° fe (""(i) ' 47 ( 2 ) ' CT ( 3 ) ' p ' q) t/ ^T( cr (i) ' er (2) > CT (3) ) p, q, p',q') 

m 

xOv(<T(j), «r (3 ), o-(i),p", q") (1 + t(2, 3) • <x ( i))|x m ) , (28) 
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G^ T ( P q,pqp'q',p"q") = 

2E^ m l° t Ki)' ff (2)> ff (3).p.q) v t ( t'k (T (i)> ,T (2).' T (3) I p I q,p' I q') 

rn 

xO k , (<r (3) , <r (2 )ff(i)p", q") (1 - <r(2, 3) • <x (1) )|x m ) . (29) 

For the relatively simple chiral NNLO 3N force [15], which we use in this 
paper, the coefficients E^J , G^'J and G' fe **' T from Eq. , (J2HJ) and (J2HD 
can be easily calculated for the whole 3N force. 

The only new coefficients, L^jk', are due to the GqV \4>) term in Eq. 

L feife '(pq J p , p",p / "q"') = 

l° fc ' °"( 2 ) ' °"(3) ' P' ^H' ( CT (2) » CT (3) > P', P") 

m 

0^(<T (1) ,<T (2) ,<T (3) ,p"',q"')|x m ). (30) 

It is possible to further simplify Eq. (|21)|) because of the scalar nature of 
the <j) functions and coefficients that one encounters there. To this aim it is 
convenient to choose our coordinate system in such a way that 

q = (o,o,g) > 



p= ( P Vl~x 2 ,0,px), (31) 
with x — p • q. 

The Cki coefficients on the left-hand side of Eq. (121)1) form a matrix with 
the determinant 

, 2916 ,,i2 2916 , n N e 

det(C) = — |pxq| 12 =— (1-z 2 ) . (32) 

It means that the C matrix might in principle cause numerical problems. 
However, the application of the inverse matrix C" 1 to the both sides of 
Eq. (|26l) and combining the C" 1 matrix elements with the L, D, D' , £7, G and 
G' coefficients leads to a scheme, which proves numerically safe. Obviously, 
the case p x q = has to be avoided in the numerical realization. 

The 3N bound state Faddeev equation is solved by iteration [T|j], us- 
ing a Lanczos type algorithm. This requires a choice of initial amplitudes 
0jy(p, q), which guarantees the Fermi character of the 3N system. Now the 
operators 02,03,06 and O7 are odd under exchange of particles 2 and 3, 
whereas Oi, O4, O5 and Og are even. Thus the scalar functions must possess 
the following symmetry properties (see Ref. [7] for the derivation): 

4t(-P- q) - (~Y +1 4t(P: q) , i = 1,4, 5, 8 , 

0$(-p, q) - RV$(P, q) , i = 2, 3, 6, 7. (33) 

Before we embark on the description of our numerical performance, we 
will show how to obtain the operator form for the full 3N bound state. This 
actually leads to yet another numerical scheme in comparison to Ref. [7] . 
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3 The Full 3N bound state 

The total 3N state 1^) is obtained from the Faddeev component as = (1 + 
P) If the isospin projected Faddeev component tptT(p,<l) = (pql^tr) = 
(pq|((t|) T\ip) is expanded (see Eq. (|22p) using the expansion functions 

0^(p,q), our task now is to find the expansion functions for Atrip, q) = 
(pq|((4) T\Pij)). We use Eq. (O but rewrite Eqs. ([i"8 ]> -([2T ]) as 

( p q I p^p^T Ip'q') = *(p' Pi (p, q))*(q' - Qi (p, q)) A S 2^ 2 S 3 . 

(pql^lT^T I P'q') = <*(p' - P 2 (p,q))%' - Q 2 (p,q))iM3 (34) 

with 



Pi(p,q) = -\p- 


3 

4 q ' 


(35) 


Qi(p,q) = p 


1 

2 q ' 


(36) 


P 2 (p,q) = + 


3 

4 q ' 


(37) 


Q2(p,q) = -p- 


1 

2 q ' 


(38) 



Consequently, we get 




+ (-)*+*' <5(p'-P 2 (p,q))5(q'-Q 2 (p,q))^ 3 F 2 s 3 } ^ T (p', q' ) 

= Y, f »'t {p? 2 Pi 3 A>T (Pi,Qi) + H t+t ' Pf 3 P 2 5 3 V' t 'T(P 2 ,Q 2 X3}3) 

t' 

On the other hand we write Atr(p, q) directly in terms of some a[ T (p,q) 
scalar coefficients 

8 

^T(p,q)=^^(p,q)O l | X m ). (40) 

i=l 

Combining (g2J, (gU) and (gUJ) we get: 

8 

^o@(p,q) O, (cr (1) ,cr (2) ,cr (3) ,p,q) \ x m ) 

i=l 

= J2 F »'T (a S 2 ^23E^'t( P 1'Qi) 4 CT (1)^(2)^(3),Pi,Qi) Ix") 
V I i=l 

+ (-)*+*' P 1 s 3 P| 3 E4^(P 2 ,Q 2 )O i ( - (1))0 - (2) , - ( 3 )) P2 ) Q 2 ) |x m ) J(41) 
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In order to get the a[^(p, q) functions we first project Eq. (|4Tj) from the 
left with (x m |Ofc(c(i),f(2) ! <''(3),P,q) and sum over m. We encounter again 
the Cki coefficients together with the new coefficients 

W k i (pq,PiQ! ) = ^(x" l |O fe (cr (1) ,cr ( 2),cr (3) ,p,q) 

JTi 

i¥ 2 i , a3 0i(o-(i),o-(a),o-(3),Pi,Qi)|x w ), (42) 
Z kl (pq,P 2 Q 2 ) = ^(x m |Ofe(T(i),T(2),^(3) I P I q) 

m 

P? S Pi 3 O i (<T (1) ,a (2) ,<T (3) ,P 2 ,Q2)\x m )- (43) 

In order to calculate Wki (pq, PiQi ) and Zki (pq, P2Q2 ) tne following iden- 
tities can be used, again following [7] 

^12 P 2 S 3 Ofe(cr (1) ,cr (2) ,<T (3) ,p,q)|x m ) = 
-\0 k {<T {2) , tr (3)> o-(i),p, q)(l + <t(23) • <x (1) )|x m > • (44) 

A S 3 P 2 S 3 Ofe(cr (1) , C r (2) , C r ( 3 ) ,p,q)|x m ) = 
-io fc (<r (3) , <r (2) ,p, q)(l - <r(23) • <r (1) )|x m ) . (45) 

Then Wki and will be calculated as 

Wki (pq,PiQi ) = - g 5Z^ m l°' £ (' T ( 1 )' cr ( 2 )'' T ( 3 )' p ' q ) 

rn 

xOi(<7 (2)) <x (3) , <t (1) , Pi, Qx ) (1 + <r(23) ■ <r { i))\x m ) , (46) 
Z fcl (pq,P 2 Q 2 ) = -- ^](x m |Ofe((T (1) ,(T (2) ,(T (3) ,p,q) 

m 

xOi(ff (3) , <7 (1)) <r (2 ),P 2 , Q2 ) (1 - <r(23) • <T ( i))|x m ) • (47) 
Using the above definitions, Eq. (jJTJ) can be written as 

8 8 

^C r fci (pq,pq)a$(P,q) = ^ F t ^ T ^ {w H (pq, P1Q1 ) 4'r(Pl> Ql ) 

4=1 t' ( = 1 

+ R t+t ' Z M (pq, P2Q2 ) 4t(^ Q2 ) )# 

The Cfci (pq, pq ) are known analytically. Thus we can easily find the 
inverse coefficients C^ 1 (pq, pq ) , which obey 

8 

} ] ^m/t (pq, pq ) Cfci (pq, pq ) — ■ (49) 

fe=i 
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Using (|49|) we finally get 

8 8 

a$\p,q) = J2 F "'T ^C m l(pq,pq)^{^(pq,P 1 Q 1 )0^(P 1 ,Q 1 ) 

V k=l 1=1 

+ (-)*+*' Z kl (pq, P 2 Q 2 ) 4^2, Q 2 ) }50) 

It is then clear that in order to calculate affi (p, q) we need three-dimensional 
interpolations. In our numerical realization we use the cubic Hermite spline 
interpolation method, described for example in Ref. |17) . 

Now we would like to pose the question of the 3N wave function normal- 
ization. The obvious starting point is 



1 = (W | V) = (V>(1 + P) I (1 + P)i/>) =3(i/)\¥) 
The full 3N bound state jif - ) can be written as 



(51) 



l*>=E Jd 3 pjd 3 q\pq) | (t^T) 

8 

E^T (P>q) °i ( (T (i) J (T (2),cr(3),P ; q) Ix" 1 ), 



where 



^(p,q) 



(p>q) + a lr (p>q> 



(52) 



(53) 



and the (p,q) functions are given by Eq. (|50D . Because of the orthogo- 
nality of the | T) and | pq) states, inserting (|S"2"j) into (fS"Tj) yields 



1 = 3 E / d3 i E (p.q) &t (p,q) 

(X™ I Oj (o-(i),<r ( 2),<T (3) ,p,q) 0^ (er (1) ,er (2) ,er (3) ,p,q) \ X n 
where the conjugate operators 0\ ("Vi), <r<2)> ^fa)) P> Q.) are 



(54) 





= Ox, 




o 2 


= o 2 , 




ot 


= -o 3 , 




1 




ol 


1 


cr 








ol 


= o 6 , 




o| 


= o 7 , 




ol 


= o 8 . 





z<T(23)-(pxq; 



r (1) x a (23)) ■ (p x q) 
r (1) x er(23)) • (p X q) 



(55) 
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If we sum both sides of (|54p over to, then we are led to 

2 = 3 E / d *p I d ^ E $t (p. q ) $ (p> q ) <5y (pq, pq ) , (56) 

tT J J i,j=l 

where 

Cij (pq,pq) = 5^(x m |0j (o- ( i),cr ( 2),o-(3),p,q) 

m 

(o-(i),o- ( a),o- (3 ),p,q) |* m ) ■ (57) 



Now all the three functions in (|56l) are scalar, so we end up with the following 
normalization condition: 



oo oo 



i = 12 tt 2 E /*^ 2 f^i 2 fdx4r (P.?.*) ftr Gw) <5y (*X58) 

where (|3"Tj) has been used. 

Using the same Cij coefficients, we can calculate the expectation value of 
the kinetic energy, (Eki n )- We start from 

(E km ) = (¥\Ho\&) = + P)\H Q \(1 + P)V) = 3 (V | #o | , (59) 

since i?o and P commute. Because the kinetic energy operator is diagonal 
for the basis states | (t|) T) and | pq) and does not depend on any angles, 
we simply get 

e oo oo 

p 2 3q 2 



(E^) = 12tt 2 Y,J2 I d PP 2 / 

■IT A A 1 " " 



tT i, 3 = l- 
1 



to 4to 



x / dx 0$ 0, q, x) 0$ (p, 9, x) Cij (x) . (60) 



Next we calculate the expectation value of the 2N potential energy, (E 2 ^ t ): 

(E™) = (V | Vu + V 13 + V 23 | *) = 3 (V | V 23 | f ) = 3 (V | V \ V ) , (61) 
where 

(f|^|f)=E/^/ d3 9E^ ( T(p,q) 

tT ^ ^ k=l 

X E /^V /dVEA'T' fp'.q') 



x (X™ I (o- (1) ,<T (2) ,CT (3) ,p,q) 

x< P q|<f4Vl^lf4'l nip'q') 



Oi(a w ,a (2) ,a (a) ,p',q') \ X m ) ■ (62) 
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Using Eqs. (gj and (JTTJ) we get 



(pq | ((t^T\V\ (t'^T')|p'q'> = 

6 

s (q - q') <5«' v tTT> (p> p' ) w j( cr (2), ct(3), p, p') • (63) 



This leads to 



I v | w) = Y, [ d3 p I ^E tir (p,q) 
£ fd 3 p'J2^ (p'.q)^ (p,p') 

T' ^ i=l 

x (x m I Ol (o- ( i),tr ( 2),<T (3) ,p,q) 
Xiy i (<7(2),CT( 3 ),p,p / )Oi (<7( 1 ),<T(2),<T(3),p / ,q) |x" 1 ) . (64) 

Again we sum both sides of Eq. over m and can write (Of \ V | <P) in 
terms of a new set of scalar coefficients, I/fc^: 

(*\v\v) = ^Y, [ d3 p f ^E (p.q) 

tT ^ ^ fc=l 

„ 6 8 

E / E (P',q) V$t- (P,P') i^(pq, P p',p'q; 

tT' 7 = 1 4=1 



(65) 



IT 

where 



L fc , i (pq ) p'p",p'V") = 

E^" 1 ! !^ 1 )' CT ( 2 )' ^W'P' ( iH( ,T (2)'' T (3)'P'' p") 

m 

xOi( - (1) , - (2))0 - ( 3 ) ,p"',q'")|x ro )- (66) 

Taking into account the scalar nature of all the functions that appear in 
Eq. (|65[) . we arrive at the final expression: 

oo oo 1 

(if I V I = 4tt 2 ^ / dpp 2 / dgg 2 / 

tJ -1 

x 53 / dp' P ' 2 / ^ /^'EEE 

T' q 10 fe=1 J =1 i=1 



Xy8 t T' (p'> 9, z') «tTT' (P>Z»'> 2/) ifcji(Pi ^p', 0' ) • (67) 
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Here we use Eq. (l3Tj) and assume that the p' vector, which appears in 
Eq. is given as 

p' = (p' Vl - x' 2 cos 4>' , p' Vl - x' 2 sin <[>' , p'x' ) , (68) 

which yields 
x' = p' • q , 



y = p • p' = m' + \/l — a; 2 vl — a/^ cos 0' . (69) 

Note that the 3N operators from Eq. (|2"B")) are given in terms of unit vectors, 
so no dependence on the magnitude of the q vector appears in the L^ji 
coefficients. 

Let us consider also the last part of the 3N Hamiltonian and calculate 
the expectation value of the 3N potential energy, (E^ t ): 

(E™ t ) = (<f | + V {2) + t/ (3) | <F) = 3 (& | V {1) | V) , (70) 



where 



tT k=l 
t'T 1 J J i=l 



X 

t' 



x(x™ I Ol (o- (1)) o- ( 2),o- (3) ,p,q; 



(pq T ' 



' 'P'q'! 



x 



Oi fai^aj.^p'.q') Ix m >- (71) 



With (USD (pq I ((4)T|yW I (t'i)T') | p'q') = ^ V#>, (p, q, p', q' ) , 
we write Eq. ([TT]) as 

(tf | yd) \$)=Yl fd 3 F f A£ $ (P'l) 
tT ^ ^ fc=l 

x E fd 3 P 'fd 3 q 'j2^Up'^') 

t> J J i=l 

x (X m I (<T(i),<T(2) ! <7( 3 ),p ! q) 

x V$ T (p, q, p', q' ) (<r (1) , <r (2) , o- (3) , p', q' ) \ X m ) ■ (72) 



14 



Repeating standard steps (summing both sides of Eq. (|72p over m and intro- 
ducing yet another set of scalar coefficients, E& T ) leads us to: 

<* i y (1) i *) = \ £ / d3 p I rf3 ?£ ^ (p.q) 

tT •* J fe=l 



X 

t 



£ / d 3 ;/ / d 3 q> 4t (p', q' ) ^' T (pq, pqp'q', p'q' ) > (73) 

V ^ i=l 



with 

^*' T (pq, pqp'q', p'q') = 

^2(x m \Ol((T(i), tT (2) , <T (3) , p, q) F t ( ^ (p, q, p', q' ) 

m 

x0. i (<T(i),<T( 2) ,<T(3),p',q')|x m ). (74) 

In the last step we use the fact that the functions in Eq. (175)) are scalar and 
reduce the number of integrals in that equation: 



OO OO 1 

(if I V {1) I !?) = 4tt 2 I d PP 2 I d 1Q 2 fdx 

j. rji \J \J \J 

n n 1 



' sm p 

* n n 



0-1 

27T OO 7T 27T 



x y dp p / d^p' sin 6 p > I d<j) p > I dq q I d9 q > sin 



7qi Sill U q i I <Mp q i 





8 

■ 

fc=l i=l 



where 

z = p' • q' = sinflp' sin# 9 < cos ((/y — <j> q >) + cos9 p i cos8 q i . (76) 

This simply follows from Eq. (pTTj) and from representing the p' and q' vectors 
in a spherical coordinate system as 

p' = (p' sin 6 p i cos (j) p i , p' sin # p < sin p / , p' cos P ' ) , 
q' = (q sin 6 q i cos g / , g' sin 9 q i sin , q' cos g / ) . (77) 



4 Numerical performance 

It is clear that our method to prepare the full 3N bound state from the Fad- 
deev amplitude, where we applied the permutation operator "to the right", 
may also be used to build iterations required for Eq. (0). One first applies 
the (1 + P) operator to the result of the previous iteration, and 
gets an auxiliary state \ip a ux) = (1 + P)\^ n ~ 1 ' 1 ) . Then the operators V 
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and are separately applied onto \ipaux) to yield \ipi) — V\^p aux ) and 
1^2) = V^\ifj aux ). The result of the new iteration, \ip^}, is just obtained as 
Go (IV'i) + IV^))- This scheme appears most efficient for the Faddeev equation 
with a 3N force and we used it in our calculations. We would like to emphasize 
that the Cij , L^ji and T coefficients introduced in the previous section 
can also be used in this new treatment of the Faddeev equation ((2J) . One can 
equally well project from the left with (x m |Ofc(f(i), f(2), f(3), P, q) or with 
(x m \Ot( (T (i)> cr (2), a (3)> P> q)- I n the latter case we arrive at 

8 

X] Ckt (pq, pq) 4t (p> q) = 

8=1 

- / ^VEEE^(P,PO%(pq,PP',p'#|T?(p'q) 

- |dVdVEE^ r (Pq.P^ 1 p'q')^ ) (p',q'). (78) 



The ^(p,q) coefficients in Eq. ((75} are the scalar coefficients for the 

(pq|((t|) T|(l + P)ip) expansion and depend on </>|^(p, q) as outlined in 
Eqs. ((50)) and ((55)) . The dependence on all relevant variables in Eq. ((T5|) can 
be made explicit as in the previous section. 

In the case where only 2N forces are included, we obtain results based on 
three different iteration schemes: 

1. |V>) = GqV{1 + P) \ip) with P applied "in both directions" as described 
in Eq. (|2BJ) of Sect. 2, 

2. \ip) = GotP \ip) as described in detail in Ref. [Jj, 

3. \ip) = G V(1 + P) as outlined above. 

Note that in all the three cases we use the same set of 3N operators (f2"5j) 
and 2N operators (0. The first set of 3N operators introduced in Ref. [7] 
combined with the 2N operators from ((5j) led to complex scalar coefficients, 
which is an unwanted feature in the bound state calculations. The numerical 
realization of all the necessary sets of scalar coefficients that appear in various 
versions of the Faddeev equation and which are necessary for the full wave 
function were prepared using software for symbolic algebra and exporting 
the expressions into the Fortran form. We used Mathematica® [18] . 

In our calculations we restrict ourselves to a chiral NNLO potential from 
Ref. [19 , with A = 550 MeV/c and A = 600 MeV/c. The operator form 
of such a potential was briefly described in Appendix C of Ref. [H] In the 
present paper we use only its neutron-proton version, neglecting the T = § 
component in the triton. 

For calculations with two nucleon forces only we use typically 36 p and 
q Gaussian points distributed in the (0,p) and (0,q) intervals, respectively. 
For the particular potential it is fully sufficient to take p « 5 fm _1 and q « 8 
fm _1 . We work with 40 x Gaussian points to represent the scalar functions 

4>tT (?) 1' x )' a tT (P> 9i x ) an d Ptr (Pi Q, x). The <f>' integrations are carried out 
using 60 Gaussian points. We will demonstrate that our calculations lead 
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to fully convergent results. As already mentioned, the Faddeev equation is 
solved by an iterative Lanczos algorithm [TB]. The initial, essentially arbi- 
trary, amplitudes fulfilled, however, the symmetry properties (|33[) . We found 
that 20 iterations are fully sufficient. A typical job with 4096 processes for 
20 iterations requires in this case less than 15 minutes on the IBM BLUE 
GENE/P computer of the JSC, Jiilich, Germany. Iterations are in fact very 
fast - one can prepare and store the essential parts of the iteration kernels 
and deal for each iteration with only two dimensional integrations. 

In the much more time-consuming calculations employing a chiral NNLO 
3N force [T5] where we deal with the integral kernel with six fold integra- 
tions we use slightly smaller numbers of p, q and x points: 32, 32 and 32, 
respectively. A corresponding job with 20 iterations requires a few hours on 
the same number of processes. In our calculations the four fold angular in- 
tegrations (see for example Eq. (fTSfl are performed using altogether rj 2500 
Gaussian points, which might leave room for improvement with respect to 
the accuracy achieved in the case without a 3N force. 

In order to save computer resources, having at our disposal the PWD 
version of the bound state code, for a given 3N Hamiltonian, we first solve 
the Faddeev equation and find the triton binding energy using a partial wave 
representation. Then for this energy we run our three dimensional calcula- 
tions and look for the eigenstate corresponding to an eigenvalue r\ closest 
to 1. Finally, for the three dimensional approaches, we construct the scalar 
coefficients for the full wave function, normalize them as in (I58|) and calcu- 
late expectation values of the kinetic and potential energies as described in 
Eqs. ([60|). (|67|) and ((751) . Results for the calculations with the chosen chiral 
NNLO 2N potential are given in Table [TJ The binding energy for this case, 
given by the partial wave based calculations, is £^=-8.299 MeV. 

Table 1 The closest to 1 eigenvalue (r\) of the kernel of the Faddeev equation for 
the given 3N energy i?3jv=-8.299 MeV, expectation values of the kinetic ({Eki n }) 
and 2N potential ({E*%)) energies in the triton as well as their sum ({Eki n ) + (Ep^ t }) 
for the guiding partial wave based results and for the methods employing the three 
dimensional formulation of the Faddeev equation (see text) . 





PWD 


Scheme 1 


Scheme 2 


Scheme 3 


V 


1.0 


0.99993 


0.99983 


0.99961 


(Ekin) 


31.963 


31.959 


31.960 


31.947 


{Epot) 


-40.258 


-40.242 


-40.243 


-40.230 


{Ekin) + {E vot ) 


-8.295 


-8.283 


-8.283 


-8.283 



Results of the calculations including additionally the NNLO chiral 3N 
force are shown in Tabled Here the calculations using standard partial wave 
representation yield £^=-8.646 MeV for the binding energy. 

Even if the new results displayed in Tables [T] and [2] show a good agreement 
with standard partial wave decomposition calculations, in the following we 
show more detailed results. We concentrate first on the case without a 3N 
force. 

It is interesting to see, if the three quite different numerical schemes 
discussed above lead to the same results. They use different sets of scalar 
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coefficients and different forms of the permutation operator, P. Addition- 
ally, the second method requires the full off-shell t-matrix, while the other 
two methods directly use the 2N potential. Last not least, in the first two 
methods the symmetry properties (|33[) were checked after the solution had 
been obtained and in the third case the same symmetry properties served to 
reduce the CPU time. In Figs. Q] [2]we show several scalar coefficients for the 
full wave function, 0$, (p, q, x), first for fixed p and a; as a function of q, then 
for fixed q and a; as a function of p and finally for fixed p and q as a function 
of a;. In the last case the symmetry properties (|33[) of the final wave function 
coefficients are clearly visible. In all the shown cases and other (not shown) 
the agreement among the three methods is very good. 

Finally, we show that it is possible to reproduce the standard partial wave 
representation of the full wave function starting from the three-dimensional 
wave function expansion (|52p. To this aim we calculate the overlap 
ipqaJM; (t|) TMj<\&mmt}, where a comprises the usual set of discrete 
quantum numbers for the 3N system, a = {I, s, j, A, 1} |14j . Here I denotes 
the orbital angular momentum in the two-nucleon subsystem, s is the total 
spin and j is the total angular momentum of this subsystem. The orbital 
angular momentum of the third nucleon is coupled with its spin (i) to give 
the total angular momentum /. Finally, j and / are coupled to yield the total 
3N angular momentum J with its projection M. Using recoupling from the 
intermediate states, where the total 3N angular momentum (L) and the total 
3N spin (S) is defined, we arrive at: 

(pqaJM; (t^J TM T \&mm t ) = 

v _ , (I s j 

6 MT , mt J2 V(2j + 1)(2/+1)(2L + 1)(25+1) < A \ I 
l,s I L S J 



Y 7T Z7T 7T Z7[ 

C(LSJ;M L ,M-M L ,M) J d9 p sin 6 P J d(f> p J d9 q sin q j 



L 

X 

M L = -L 



2tt tt 2tt 

l<t>q 







:Va iM MP,?)I>$(p,q) 



i=l 

,( 1\ „ 

X 



2 

where 



s\)SM-M L \O i {a {1) ,a [2) ,a^ ) ,p,ci) \ x m ) , (79) 



3f A Mi (p, q) = C ( lXL ^ m i' M L ~ m u M L ) Y lmi (p ) Y XML - mi {q ) (80) 

m;=— ! 

C{j\jij\m\mirn) are the Clebsch-Gordan coefficients and Yi mi (p) denote 
the spherical harmonics. The essential parts of the overlap ([79")) . 

SM - M L \Oi{cT {1) ,a m ,CT m ,^q) \ x m ) , (81) 
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can again be prepared with help of the symbolic algebra software. It is also 
clear that the overlap ([?§]) can be used for further verification of the three 
dimensional solution. In particular overlaps with | pqaJM; (t^J TMt) car- 
rying quantum numbers different than the ones of 3 H should be numerically 
zero. We verified very many cases but can show only a few examples (see 
Figs. HHSJ), which demonstrate that agreement with standard partial wave 
decomposition is perfect not only for the dominant partial waves but also 
for the smaller triton components. For this comparison we used the three 
dimensional wave function obtained with the third method. 

Although the numbers of p and q points in our calculations including 
the chiral NNLO 3N force were smaller than in the case where only 2N 
force were present and the four fold angular integration in the integral kernel 
including the 3N force might be improved, we show in Figs. [5] and [7] that the 
projections of the full wave function onto the partial wave states agree quite 
well with the standard partial wave representation of the 3N bound state. 
This is remarkable, since the partial wave calculations were based on Eq. (fTJ) 
and the three dimensional calculations used Eq. ©. 

It is important to notice that the normalization condition (analogous to 
Eq. dSHD) 

g OO OO 1 

1 = 47f2 £ £ [ d PP 2 [ d 11 2 [ dx $T (P.9>*) (P. 9.*) &ii (*)(92) 
leads to the following normalization of the partial wave projections 

oo oo 

Y, Jdpp 2 Jdq q 2 ((pqa \ <P) f = 1 , (83) 
a 

and provides an additional check of numerics. 



5 Summary and Outlook 

Using an expansion of the Faddeev amplitudes in products of scalar functions 
04t(P' °0 with products of scalar expressions in spin and momenta as well as 
corresponding expansions of the 2N potential, 2N i-operator and 3N forces, 
the Faddeev equations for the 3N bound state can be reformulated into a 

strictly finite set of coupled equations for the amplitudes (f>\j,(p, q), which 
depend only on three variables p , q and p ■ q. 

In this paper we discuss various forms of the Faddeev equation in the 
operator form. In comparison to Ref. [Jj we propose two new schemes which 
use directly a 2N potential (without a 3N force) and a new approach to 
the Faddeev equation with inclusion of a 3N force. For each case we derive 
corresponding sets of scalar coefficients. For the three schemes without a 3N 
force and one scheme including a 3N force we provide numerical realizations 
based on the chiral NNLO 2N [T9lfT2] and 3N [15] forces. 
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Table 2 The closest to 1 eigenvalue (rf) of the kernel of the Faddeev equation for 
the given 3N energy i?3jv=-8.646 MeV, expectation values of the kinetic ({E^in}), 
2N potential 3N potential ({E™ t }) energies in the triton as well as their 

sum ((Ekin)+(E$%)+(E$%)) for the guiding partial wave based calculations and 
for the method using the three dime nsio nal formulation of the Faddeev equation 
with a 3N force (see text around Eq. l|78|lV 



PWD 



~3D~ 





1.0 


0.99976 


(Ekin) 


33.448 


33.412 


(Ep t) 


-41.329 


-41.273 


(Epot) 


-0.765 


-0.770 


{Ekin) + {E pot ) + {E pot ) 


-8.646 


-8.631 






0.005 r 




- 




-0.005 - 




-0.01 - 




-0.015 - 




-0.02 - 


d" 

d. 


-0.025 - 




-0.03 - 




-0.035 - 




-0.04 7 




-0.045 ' 




-0.05 L 




0.005 




|[fm ; 


-0.005 - 


X 

d- 


-0.01 - 


d. 






-0.015 - 




-0.02 7 




-0.025 - 




2 3 4 
q [1/fm] 

i=5,t=1 






0.0002 r- 




- 


"e 


-0.0002 | 




-0.0004 I 




-0.0006 - 


d. 


-0.0008 - 




-0.001 - 




-0.0012 - 




-0.0014 - 




Fig. 1 Selected scalar coefficients /3 tT ! (p, q, x) for T = 

function in the three dimensional form for p w 1.4 fm 
of q obtained with the first (dotted line) , second (dashed line) and third (solid line) 
iteration scheme without a 3N force (see text). All three lines essentially overlap. 



| for the full 3N wave 
1 and x w 0.55 as a function 
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The present paper is a continuation of Ref. [7] and presents further el- 
ements of the three dimensional formalism. In particular we show how to 
obtain the total 3N bound state from the Faddeev component, using another 
form of the permutation operator. The paper provides also formulas for the 
bound state normalization and the expectation values of all the parts of the 
3N Hamiltonian. 

We show our numerical results for the binding energies and display ex- 
amples of the scalar expansion functions for the full 3N bound state. We 
demonstrate stability and reliability of our three dimensional treatment of 
the bound state Faddeev equation by comparing results obtained with the 
three different methods. One of the methods required a new fast way to ob- 
tain the full off-shell i-matrix in the operator representation. We will report 
on these calculations in a separate paper. We give all relevant information 
about the number and distribution of the p, q and x points used in our cal- 
culations. We are ready for benchmark calculations of the 3N bound state in 
the three dimensional formalism. Our final results for the full wave function 
obtained in the three dimensional schemes compare very well with the stan- 
dard results based on the partial wave decomposition. This is the best proof 
of the stability and accuracy of our numerical performance. 

The next step is the extension of our calculations to include the Coulomb 
force in 3 He. The inclusion of this long-ranged force might be easier in the 
three dimensional formalism. We also plan to apply a three-dimensional oper- 
ator approach to the problem of 3N scattering and use the formalism outlined 
in [2D]. 
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Fig. 2 Selected scalar coefficients /3 

function in the three dimensional form for q 
function of p. Lines as in Fig. [T] 



for T = i for the full 3N wave 



0.72 fm" 1 and x 



-0.61 as a 
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i=3,t=1 





0.0004 r 




0.00035 




0.0003 




0.00025 


dr 
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Fig. 3 Selected scalar coefficients {3$ (p, q, x) 

function in the three dimensional form for p ; 
function of x. Lines as in Fig. [T] 



for T = 
i 2 fm~ and q 



| for the full 3N wave 
2.6 fm _1 as a 
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Fig. 4 The partial wave projected full 3N wave function (pqa^; (t^) ^l^l \) 
for p « 0.44 fm _1 and four different partial wave states a: a = 1 corresponds to 
I = 0,s = 0,j = 0, A = 0,1 = \,t = 1, a = 5 gives i = 0, s = 1, j = 1, A = 0, / = 
i, t = 0, a = 18 comprises i = 2, s = l,j=2, A = 2, / =|,t = and for a = 26 we 
have I — 3, s = 1, j = 3, A = 3, / = |, t = 1 as a function of <?. The results obtained 
directly using the par tial wave representation (dashed line) are compared with the 
results based on Eq. (|79|l . Two lines practically overlap. 




0.5 1 1.5 

p[1/fm] 




Fig. 5 The partial wave projected full 3N wave function (pqa~; ^||^'^^) 
for q « 0.9 fm~ and the same four different partial wave states a as in Fig. Uas 
a function of p. 
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a- 1 a- 5 




p[1/fm] p[1/fm] 

Fig. 7 The same as in Fig. [5] but with inclusion of the chiral NNLO 3N force in 
the both types of calculations. 



